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ABSTRACT 



Context. Optical observations from 1989 of the Small Magellanic Cloud (SMC) B[e] supergiant star LHA 115-S 23 (in short: S23) 
CIh' revealed the presence of photospheric Hei absorption lines, classifying S 23 as a B8 supergiant. In our high-resolution optical spectra 

from 2000, however, we could not identify any Hei line. Instead, the spectral appearance of S 23 is more consistent with the classifi- 
cation as an A 1 supergiant, maintaining the so-called B[e] phenomenon. 

Aims. The observed changes in spectral behaviour of S 23 lead to different spectral classifications at different observing epochs. The 



o 



aim of this research is, therefore, to find and discuss possible scenarios that might cause a disappearance of the photospheric Hei 
^ i ' absorption lines within a period of only 1 1 years. 

Methods. From our high-resolution optical spectra, we perform a detailed investigation of the different spectral appearances of S 23 
based on modern and revised classification schemes. In particular, we derive the contributions caused by the interstellar as well as 
the circumstellar extinction self-consistently. The latter is due to a partly optically thick wind. We further determine the projected 
qq ' rotational velocities of S 23 in the two epochs of spectroscopic observations. 

Results. Based on its spectral appearance in 2000, we classify S 23 as Al lb star with an effective temperature of about 9000 K. This 
, classification is supported by the additional analysis of the photometric UBV data. An interstellar extinction value of E(B — V)— 0.03 

■ is derived. This is considerably lower than the previously published value, which means that, if the circumstellar extinction due to 

the stellar wind is neglected, the interstellar extinction, and hence the luminosity of the star, are overestimated. We further derive a 
rotation velocity of v sin i ^ 150 km s , which means that S 23 is rotating with about 75% of its critical speed. The object S 23 is thus 
the fourth B[e] supergiant with confirmed high projected rotational velocity. The most striking result is the apparent cooling of S 23 
by more than 1500 K with a simultaneous increase of its rotation speed by about 35% within only 1 1 years. Since such a behaviour is 
excluded by stellar evolution theories, we discuss possible scenarios for the observed peculiar variations in S 23. 
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5r i Introduction metric observations, e.g., by Magalhaes (1992), Magalhaes et al. 



3), and Melgarejo et al. (2001 ), in rare cases even from op- 

B[e] stars belong to one of the puzzling phenomena in modern tical imaging with the Hubble Space Telescope, as wa s the case 
astrophysics. While many of their members have been known for for the g alactic B[e] star Hen 2-90 (Sahai et al. 120021 Kraus et 
decades, the number of new detections of B[e] stars is growing a i. [2005), and from interferometric observations with, e.g, the 
tremendously (see, e.g., Kraus & M iroshnichenko 2006). The VLTI/AMBER and VLTI/MIDI instruments (e.g., Domiciano 
classification by Lamers et al. d 19981 ) of B[e] stars according to de Souza et al. 12007] ). An observed strong infrared (IR) excess 
their evolutionary phase made clear that B[e] stars occur in the emis sion due to warm and ho t circumstellar dust (e.g., Z ickgraf 
pre- as well as in the post-main sequence phase, and the most e t al. [19851 [19861 [19891 [19921 Gummersbach et al.QM5j Borges 
popular class is formed by the B[e] supergiants with confirmed Fernandes et al. 120071 ) indicates the presence of a high-density 
members predominantl y in th e Magellanic Clouds (see the re- disk. This disk must furthe r be th e location of the neutral ma- 
cent review by Zickgraf [2006]). terial like Oi (Kraus et al. [20071 ) and the molecular gas, ex- 
The derivation of stellar parameters of B[e] stars is com- pressed by C O first-ove rtone b ands (Mc Gregor et al. 1988a, 
plicated because in most cases the stellar photospheres are hid- 1988b! [T9%9l Morris et a l.|1996l Kraus et al. |2000| ) and TiO band 
den inside the dense and optically thick winds. Thus, the spec- emission (Zickgraf et al. |19891 >. The proposed scenario for B[e] 
tra are contaminated with wind continuum emission in the form supergiants thus suggests that their circumstellar material con- 
of free-free and free-bound emission as well as with an enor- sists of a disk-like structure formed by a high density equatorial 
mous amount of emission lines that are generated in the cir- wind with low outflow velocity, and a normal B-type line-driven 
cumstellar wind material, which is usually non-spherically sym- wind of high v elocity and muc h lowe r density in polar directions 
metric. Evidence for non-sphericity thereby comes from polari- (Zickgraf et al. 1985, Zickgraf 2006). A possible mechanism for 
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the formation of outflowing high density disks might be rapid ro- 
tation of the central star (see, e.g., Bjorkman & Cassinelli 1993; 
Cure 2004; Cure et al. 2005). In addition, if the star is rotating 
rapidly, the wind material in the equatorial plane recombines in 
hydrogen in the vicinity of the stellar surface, as has recently 
been shown by Kraus (2006 ). And for at least three B[e] super- 
giants projected rotational velocities at a significant fraction of 
their critical speed have indeed been found (Gummersbach et al. 
119951 Zickgraf l2000t 120061 ). hinting to the possible connection 
between the presence of a circumstellar disk and stellar rotation. 

For the Small Magellanic Cloud (SMC) star LHA 1 15-S 23 
(AzV 172, 2MASS J00555380-7208596, in the following ref- 
ered to as S 23) that is investigated in this paper, no rotational 
velocity has been derived yet. In addition, not much is known 
about the orientation of the system, although it has been sug- 
gested to be seen more or less edge-on, based on the low ex- 
pansion velocity derived from the Balmer lines (Zickgraf et al. 
1992 referred to as ZSW in the following). This object was 
first classified by Azzopardi & Vigneau (1982) as a B8 super- 
giant, and the detection of several photospheric Hei absorp- 
tion lines in the optical mid-resolution spectrograms of ZSW 
is in agreement with the classification of S 23 as a late B-type 
supergiant. These authors further fixed the interstellar extinc- 
tion at E(B - V) = 0.10, derived an effective temperature of 
T e s = 11 000 K, and a surface gravity of log g - 2.0. This led to 
M bo i = -6.4, log L«/L = 4.46,/?* = 45 R Q , and the assignment 
of the luminosity class lb. The spectra of ZSW showed further 
strong Balmer line emission, [Oi] emission, as well as many Fen 
and [Fen] emission lines. In combination with the pronounced 
infrared excess in the photometric data, which can be ascribed to 
warm or hot circumstellar dust, ZSW assigned S 23 to the group 
of B[e] supergiants. 

About eleven years later, during our observing mission of 
Magellanic Cloud B[e] supergiants, we observed S23 using a 
high-resolution optical spectrograph. Our spectra confirm the 
presence of the strong Balmer emission lines as well as of the 
[Oi], [Fen], and Fen emission lines. But interestingly, even with a 
more than 10 times higher resolution, our spectra do not display 
any Hei line either in absorption nor in emission. The absence of 
the Hei lines speaks more in favour of a cooler central object than 
a late B supergiant. This paper is therefore aimed to investigate 
in detail the spectral appearances of S 23 at the two epochs of 
observations and to discuss possible scenarios that might cause 
the observed peculiar variations of S 23. 

The paper is structured as follows: our observations are 
described in Sect.|2] Then, we perform the spectral classifica- 
tion of S 23 for the two different observing epochs, using mod- 
ern, improved classification schemes (Sect.O and determine the 
contributions of both, interstellar and circumstellar extinctions 
(Sect.|4j>. From the detected photospheric Mgn absorption lines, 
we next derive the projected rotation velocity of S 23 at the two 
epochs of observations (Sect.|5j. In Sect.|6] we discuss possible 
scenarios that might explain the observed peculiarities, before 
we summarise our major results in Sect. [7] 

2. Observations 

We obtained the optical spectra with the 1.52-m telescope at the 
European Southern Observatory in La Silla (Chile) using the 
high resolution Fiber-fed Extended Range Optical Spectrograph 
(FEROS). FEROS is a bench-mounted Echelle spectrograph 
with fibers, which covers a sky area of 2" of diameter, with a 
wavelength coverage from 3600 A to 9200 A and a spectral res- 
olution of R = 55 000 (in the region around 6000 A). We adopted 
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Fig. 1. Parts of the rectified FEROS spectrum (dashed lines) 
showing the regions around several expected Hei lines. For bet- 
ter visualization, we convolved the FEROS data to a resolution 
of R = 4900 as was used by ZSW (solid lines). 



its complete automatic on-line reduction. The high-resolution 
spectrum of S 23 was taken on October 14, 2000, with an ex- 
posure time of 3600 seconds. The S/N ratio in the 5500 A region 
is approximately 40. Since within the same night we were able 
to take two consecutive spectra of the star, which do not show 
significant differences, we added them up for a better S/N ratio 
achievment. 

In Fig.[T] we display parts of the rectified high-resolution 
FEROS spectra. The first three panels show the regions around 
the positions of several Hei lines, which were identified by ZSW 
in their 1989 low-resolution (i.e. R = 4900) observations (see 
Fig. 4 of ZSW), with the Hei /14472 line being the strongest one 
with almost equal equivalent width to the adjacent Mgn /14481 
line. This Hei line is clearly absent in our FEROS spectrum (or 
at least below the detection limit) as well as the other two Hei 
lines at /1/14388, 4438, whose detection was previously reported. 
The fourth panel of Fig.[T] displays the region around the Hei 
/15876 line. If Hei is present in a stellar spectrum, then this line 
is usually quite strong. But obviously, it is also absent from our 
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Fig. 2. Balmer lines in the spectrum of S 23. 



Table 1. Mean heliocentric radial velocities derived from differ- 
ent ions compared to those derived by ZSW. 
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FEROS spectrum. We checked the complete FEROS spectrum, 
but we could not identify a single Hei line. 

Due to our rather low S/N ratio and the fact that the pho- 
tospheric Hei absorption lines are expected to be broad and 
shallow, the high-resolution spectrum might not be able to dis- 
play these lines properly. Instead, for these lines a (much) lower 
spectral resolution might be advantageous as has been shown 
by Verschueren d 1 99 li t. We therefore convolved our data to 
the same resolution as the ZSW spectrum and overplotted the 
low-resolution spectrum in Fig.Q] However, even with the low- 
resolution data, the Hei lines do not show up. 

While other known B[e] supergiants do not display even a 
single photospheric absorption line in their optical spectra due 
to their high density and therefore optically thick winds, S 23 
displays at least one clearly detectable line, namely the Mgn 
doublet at /14481. And also the higher Balmer lines (HS, Hy, 
see Fig.|2j, as already mentioned by ZSW, show indications for 
a photospheric absorption component on which PCygni pro- 
files seem to be superimposed. Some other lines in our spectra 
also display absorption components, like Fen, Tin, and Sin (see 
Fig.Q]). These absorption components are, however, much nar- 
rower than the Mgn line, indicating that they are circumstellar in 
origin. 



We measure the heliocentric radial velocity for several lines 
in our spectra: from the narrow pure emission lines of the for- 
bidden transitions (like [Oi], [Fen], see Fig.Q]), from the narrow 
absorption components observed for most of the permitted tran- 
sitions (like Fen, and two lines of Sin and Cm), and from the 
narrow absorption components of the Balmer linefl For each 
element we derive mean values, which are listed in TableQ] For 
comparison, we also included in this table the values obtained 
by ZSW. While our metal lines all deliver values in the range 
155 . . . 160 kms -1 , the values of ZSW indicate a much larger 
scatter between the radial velocities of the forbidden and the 
pemitted lines, which had been interpreted by these authors as 
a slow outflow component. In fact, the narrow absorption com- 
ponents of the permitted lines, often overlaid on broader emis- 
sion components, are not of pure atmospheric origin. Instead, 
they are generated due to absorption by the circumstellar mate- 
rial. The same holds also for the narrow absorption components 
seen in the Balmer lines. For a proper radial velocity determina- 
tion we therefore rely on forbidden lines, which display narrow 
but symmetrical line profiles. From these lines, we thus obtain 
a mean radial velocity of v m d = 156.5 ± l.Okms -1 , which we 
will use here. It is slightly lower than the value found by ZSW 
from their [Oi] lines, but this might be a spectral resolution effect 
rather than a stellar effect. 

The largest difference between the ZSW and our data is 
reached by the two emission peaks of the Ha line. Here, a clear 
shift in velocities is present. This shift is about 40 kms-' for the 
blue peak, while it is only 10-15 km s _I for the red peak. We 
tested the influence of the different resolutions by smearing our 
Ha line to the same resolution as ZSW. The measured radial ve- 
locities are in much better agreement with the ZSW values, so 
that we can conclude that the differences in the appearance of the 
Ha line in ZSW and in our spectrum are most probably caused 
by the different spectral resolutions. 

3. Spectral classification of S23 

3.1. Temperature classification 

The determination of the spectral type and luminosity class of 
stars with metallicities different from the galactic one has been 
found to be a difficult task in the past. The reason for this is 
the fact that with decreasing metallicity the metal lines in the 
spectra become very weak and hard to detect. The proper de- 
termination of the spectral type is rather complicated, especially 
for B-type stars which do not or only very weakly show Hen 
lines in their spectra. For these objects, the classification has to 
be performed with the help of the metal lines and is, for galac- 
tic objects, normally based on the line strengths of silicon and 
magnesium, usually in comparison to helium. As shown, e.g., 
by Lennon et al. Q993 ) for galactic B supergiants, the equiva- 
lent widths of many spectral lines show a typical behaviour with 
spectral type, i.e. effective temperature. And later on, Lennon 
( I1997D suggested determining the spectral classes of B-type su- 
pergiants in the SMC according to line ratios in such a way that 
the same trend in line ratios with spectral type is achieved as for 
solar metallicity (Lennon et al. |1992| ). His classification scheme 
(see Table 1 of Lennon, 1 19971 1, which has been extended to G- 
type stars by Evans et al. ([2004), is used by many authors to 



1 For the Balmer lines, we list the mean radial velocity derived from 
the central absorption component for 4 lines (H/3, Hy, H<5 and He - H/3 
and He were not observed by ZSW), while for Ho- the radial velocity 
derived for the central absorption ( ca ), the red peak ( rp ), and the blue 
peak ( bp ) are given. 
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Fig. 3. Effective temperature dependence of the line ratio Mgn 
/14481 over Hei /14472 calculated from model stellar atmo- 
spheres with logg = 2.5 at solar abundances. The behaviour be- 
tween spectral type and effective temperature taken from Evans 
& Howarth (120031 1 is indicated in the bottom label for solar abun- 
dance, and in the top label for SMC abundance. The two boxes 
indicate the observed ratios with their errors for the ZSW spec- 
trum taken in 1989 and our spectrum taken in 2000, respectively. 

classify B-type supergiants in many galaxies with abundances 
similar to the SMC value (e.g., Evans & Howarth 2003 ; Bresolin 
et al. 120061 Bresolin et al. l20l)7i l. 

The spectral type of S 23 has been determined by Azzopardi 
& Vigneau d 19821 1 as B8, and later on, ZSW, derived an effec- 
tive temperature of about 1 1 000 K based on comparison of the 
observed spectral energy distribution to Kurucz d 19791 1 model 
atmospheres. Inspection of their mid-resolution spectrum indi- 
cates the presence of the Hei AAA12 and Mgn /1448 1 lines, with 
the Mgn line being slightly stronger than the Hei line. Lennon et 
al. ([1993 ) already mentioned the sharp decrease in Hei equiva- 
lent width from early to late B-type supergiants, and the simulta- 
neous strong increase in Mgn /1448 1 equivalent width. The rela- 
tive strength of these two lines is therefore one of the main clas- 
sification characteristics of especially late B-type supergiants 
(e.g., Bresolin et al. 2006J Lennon[T997|>. According to the clas- 
sification scheme of Lennon (119971 ), the line ratio as inferred 
from the ZSW spectrum would thus suggest a spectral type be- 
tween B8 and B9 of S 23 in the year 1989. 

Our own high-resolution spectra display only the Mgn /1448 1 
line while the adjacent Hei /14472 line is definitely below the de- 
tection limit. This would speak more in favour of a star with 
spectral type AO or later. To quantify these classifications and 
to constrain the possible range in effective temperature at the 
two observing epochs, it is necessary to calculate theoretical 
line ratios of Mgn /1448 1 over Hei AAA12. This is done by us- 
ing the code SYNSPEC (see Hubeny & La nz 120001 ) to compute 
line identification tables based on Kurucz d 19791 ) model atmo- 
spheres in local thermodynamical equilibrium (LTE) and Kurucz 
J 19931 ) line lists. These calculations have been performed for ef- 
fective temperatures ranging from 8 000K up to 15 000K and 
solar abundances. With these line identification tables that con- 
tain information about the expected line equivalent widths, we 
then calculated the Mgn /I4481 over Hei /14472 line ratios. The 
results are shown as the dashed line in Fig. [3] To convert the ef- 



fective temperature scale into a spectral type, we made use of 
the recent classification of Evans & Howarth (2003), who re- 
defined effective temperatures for galactic (indicated along the 
bottom axis of Fig. [3]) and SMC supergiants (indicated along the 
top axis). In addition, we follow the approach of Lennon ( 119971 ), 
who suggested that the same line ratio mirrors the same spectral 
type. We find that the transition from a line ratio smaller than 
one to a ratio larger than one happens at an effective temperature 
of about 12 000 K for solar abundances, which corresponds to a 
spectral type between B8 and B9, in agreement with what has 
been claimed by Lennon ( 11997I ). 

Concerning the measured line ratios of S 23 we could derive 
a value of 1 .25 ±0.25 from the low-resolution spectrum of ZS W3 
This value, with its error, is included in Fig. [3] Its location in the 
diagram indicates that in 1989 the spectrum of S 23 was that of a 
B8.5-B9 star with effective temperature T eff =* (1 1 000 ± 500) K. 
This classification is in fair agreement with the one derived by 
ZSW. 

From our own spectrum taken in 2000, we derive a line ratio 
of 2.85 ± 0.25. This ratio has to be considered as a lower limit 
since we used, for the equivalent width of the Hei AAA12 line, 
the value we consider the detection limit. With this line ratio, we 
therefore determine an upper limit to the effective temperature 
that is mirrored by this line ratio. We find, for S 23, a spectral 
type later than AO and an effective temperature % 9 500 K. 
If the spectral appearance of S 23 at the two observing epochs 
is indeed determined by the temperature of the central star, then 
the change in line ratio suggests that the star must have cooled by 
at least 1500 K within a period of 1 1 years only, a scenario that 
is definitely excluded by stellar evolution theory (see Sect. 16. It . 
In the following, we will therefore refer to this behaviour as an 
apparent cooling of the star. 

3.2. Luminosity class 

For the SMC, great effort has been undertaken for a proper lumi- 
nosity classification of the supergiant population. In principle, 
several methods for the determination of the stellar luminosity 
exist and have been successfully applied in the literature. For 
O-type to early B-type supergiants, the Hei/Hen line ratios are 
good and, especially, metallicity-independent luminosity indica- 
tors (e.g., Walborn 1983 ). For B and later type supergiants, the 
Hy equivalent width is an ideal and metallicity-independent tool 
for the luminosity class determination, as shown by Azzopardi 
(119871 ), while the use of metal lines as suggested by Massey et 
al. ([T995 ) delivers unreliable results when applied to stars in 
low-metal content galaxies like the SMC (see the discussion in 
Lennon [T997l 

For S 23, the luminosity class determination is complicated 
due to fact that its Hy line shows an additional emission com- 
ponent (see Fig. |2]). In this case, the equivalent width method of 
Azzopardi ( 1987 ) is not applicable. 

ZSW had derived the luminosity class lb from determining 
E(B - V) based on the photometric B and V data of S 23 from 
Azzopardi & Vigneau (1982). Using the intrinsic colors for a 
B8 supergiant from Schmidt-Kaler (1982), they found a value 
of E(B — V) = 0.10. However, as we will show in Sect.[4] this 
extinction determination is based on assumptions that are not re- 
liable for typical B[e] stars, and consequently, the E(B - V) value 



2 This value has been derived by hand from the published spectrum 
because the digital spectra are not available anymore (Zickgraf, private 
communication). 
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Table 2. Photometric data of S 23 (in mag) collected from the literature and sorted by their year(s) of observation. 
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Year 
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V 


R 


J 


H 


K/K s 


References 


1972-74 


13.27 


13.43 


13.37 










Azzopardi & Vigneau ( 1982) 


1988 










13.04 


12.79 


12.16 


Zickgraf et al. 09921 


1998 










13.02 


12.76 


12.15 


Skrutskie et al. (120061 2MASS) 


1999 


13.11 


13.38 


13.31 


13.19 








Massey (2002) 



derived by ZS W overestimates the real amount of interstellar ex- 
tinction and, therefore, the luminosity of S 23. 

From our FEROS spectrum alone, we are not able to make 
any estimates about the stellar luminosity. We therefore searched 
the literature for complementary photometric data. These data, 
which have usually errors on the order of 0.01 mag, are sum- 
marised in Table|2l and the K s band belongs to the 2MASS ob- 
servations. Inspection of this table delivers the following impor- 
tant information: 

- The UBV data of S 23 have changed substantially between 
the observations taken in 1972-74 and those taken in 1999, 
while the near-IR photometry appears constant in the mea- 
surements taken in 1988 and in 1998. 

- The UBV data of Azzopardi & Vigneau (|1982| l were taken 
14-16 years before ZSW took their near-IR data. With the 
knowledge of a possible variation/cooling of the central star 
as discussed above, the combination of these two datasets for 
the construction of the spectral energy distribution is there- 
fore most certainly unreliable. 

- In 1999, the star appears to be more luminous in the UBV 
bands than it was during the 1972-74 observing period. 

- The UBV data of Massey (2002) were taken only one year 
earlier than our FEROS spectra, and we use them as comple- 
mentary data to our optical spectra for a tentative luminosity 
determination. 

From the photometry of Massey (2002 ) we can estimate the 
luminosity class of S 23 in the year 1999, if we completely ne- 
glect the interstellar extinction towards S 23 for the moment. 
From the distance modulus of /j = 18.93 (Keller & Wood 
2006), we thus find a first-guess absolute visual magnitude of 
My = -5.62, assigning to S 23 the luminosity class lb. 

4. Interstellar and circumstellar extinction 

In the previous section, we determined the luminosity class of 
S 23 as lb. For this, we ignored the presence of any interstellar 
extinction value, knowing that this delivers only a lower limit 
to the real stellar luminosity. In this section, we therefore intend 
to derive the values of the interstellar and possible circumstel- 
lar extinction. This is necessary for a proper dereddening of the 
photometric data and a reliable stellar luminosity determination. 

To derive E(B - V) from the color index of a star, one should 
not only rely on the B — V color index, but also check, if pos- 
sible, the result by using a different color index, e.g. U - B. 
The observed colors from Azzopardi & Vigneau (1982) given 
in Table|2] serve as an example. According to the classification 
of S23 as B8Ib, ZSW derived a value of E(B - V) = 0.10 
from (B - V). Using as well the color index (U - B), and 
the intrinsic color of Johnson ( 1 19661 ), we find an extinction of 
E(B — V) = 0.49. For the calculation of this value, we used the 
relation E(B - V) — E(U - B)/0.11, which has been found from 
observations by Leitherer & Wolf ( 1984). The value found from 
the U - B color index thus seems to indicate an extinction five 
times higher. 



One might argue that, since the star seems to be somehow 
variable, the spectral type derived from the optical data of ZSW 
cannot be applied to the UBV data taken about 14 years earlier. 
Therefore, we checked the possible interstellar extinction val- 
ues derived from the more recent set of UBV data from Massey 
(2002), using the spectral classification derived from our optical 
spectra that have a time offset of about one year, only. 

From Fig. [3] we found the spectral type of S 23 as AO or 
later, and in the following we will deal with two stellar scenar- 
ios, namely an AO lb star with T e s = 9500 K, and an A 1 lb star 
with T e ff — 9000 K. Using the tables with intrinsic colors from 
Johnson ( 1966 ), we find the following extinction values: 

E(B - V) = | [^g J from {(^Ig)} for AO lb, 
and 

= from {$-¥>} for A1Ib - 

Again, both methods (even considering an error in the UBV band 
fluxes of about 0.01 mag) deliver results that differ by a factor 
2-3. This discrepancy in extinction implies that, obviously, the 
interstellar reddening is not the only extinction source acting on 
the observed UBV data. Instead, circumstellar material along the 
line of sight is causing an extinction of the stellar flux as well. 

Possible circumstellar extinction sources might be, on the 
one hand, circumstellar dust with a wavelength dependent ex- 
tinction behaviour different from the the global interstellar one. 
On the other hand, since B[e] stars are known to have high mass 
loss rates, the circumstellar extinction might be caused by the 
high density wind that drops partly optically thick at optical 
and/or UV wavelengths. 

To test which circumstellar extinction source might be act- 
ing on the UBV fluxes of S 23, we performed some simple test 
calculations for both scenarios. For the circumstellar dust sce- 
nario, we adopted a grain size distribution according to Mathis 
et al. d 19771 ) for either silicates or amorphous carbon grains with 
SMC abundances of Si and C taken from Venn (199 9) . The opti- 
cal depths of the circumstellar grains, resulting in a wavelength 
dependent dust extinction, are then calculated for the scenarios 
of an AO lb and an A 1 lb star. Interestingly, we found that neither 
the silicates nor the carbon grains (nor any combination of both) 
can account for the circumstellar extinction of S 23. We therefore 
conclude that circumstellar dust is most probably not the (major) 
source of the extra extinction. The location of the circumstellar 
dust around S 23, whose presence might be inferred from the IR 
photometry, should thus be out of the line of sight. 

The situation is more complicated when the extinction is 
caused by the stellar wind. In this case, the wind not only ab- 
sorbs and scatters the stellar light, instead it contributes with ad- 
ditional emission even at optical and UV wavelengths, as has 
recently been shown by Kraus et al. d2008bl ) for classical OB 
supergiants, and by Kraus et al. ( 120071 ) for the LMC B[e] super- 
giant R 126. 
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Fig. 4. Fitting process of the Mgn /i448 1 line with a pure rotation profile (left panel), a pure Gaussian profile (mid panel), and our 
best fit model (right panel) with a profile consisting mainly of stellar rotation plus some macro-turbulence. 



We test the influence of the wind to the UBV band fluxes by 
assuming for simplicity (and since nothing is known about the 
real structure and geometry of the wind around S 23) a spher- 
ically symmetric isothermal wind, and calculate the wind con- 
tribution with the model of Kraus et al. (2008b). For the termi- 
nal wind velocity we use the escape velocity, which was found 
from observations of late B-type and early A-type supergiants 
by Lamers et al. (Q/995) to be a reasonable approximation. The 
only free parameter in the wind calculations is thus the mass loss 
rate, M, which is varied over a certain range to calculate the wind 
extinction and emission in the UBV bands. 

Interestingly, we could not constrain the circumstellar ex- 
tinction for the scenario of an AO lb star, while the scenario of 
an A 1 lb star delivered reasonable values for both the mass loss 
rate, with M =* 3 x 10~ 6 M o yr _1 , and the interstellar extinction, 
with E(B — V) — 0.03, constraining the classification of S 23 
more towards an Al star rather than an AO star. 

Compared to the value found from the observed B — V color 
index, our value is slightly lower, indicating that the neglect of 
the wind influence to the emergent flux from the atmosphere 
leads to an overestimation of the real interstellar extinction. Such 
behaviour has also recently been found for classical B-type su- 
pergiants (Kraus et al. 2008a). 

With the tentative classification of S 23 as Al lb star with an 
effective temperature of roughly 9000 K, and with the derived 
interstellar extinction of Ay 0.10, we find the following set 
of stellar parameters: log(L„/L ) 4.3, 7?» ^ 55 — 60 R Q , and, 
from the comparison with stellar evolutionary tracks calculated 
for SMC stars by Charbonnel et al. (|T993), an initial mass of 
M im = 9.5...11M . 

Finally, the small interstellar extinction value confirms the 
luminosity class adjustment of lb, meaning that S 23 belongs to 
the group of low-luminosity B[e] supergiants, which has mem- 
bers in the LMC (Gummersbach et al. 1995) as well as in the 
Milky Way (Borges Fernandes et al. 2007), and which should 
probably be extended at least to early A-type stars. 



5. Stellar rotation 

With the classification of S 23 as an early A-type supergiant at 
the time of our observation, we next intend to determine its stel- 
lar rotation speed given by the parameter v sin i. For late B-type 
and early A-type stars and supergiants, the most reliable pho- 
tospheric line for the stellar rotation determination is the Mgn 
A4481 doublet (e.g., Rover [20051 Royer et al. |2002l , which is 
also the only clearly detectable photospheric line in our spectra. 

It has recently been shown that the Fourier analysis provides 
an excellent tool for the rotation velocity determination even 
for hot stars (see Royer 120051 Simon-Diaz & Herrero 2007). 
However, for the Fourier analysis to work properly, it is nec- 
essary to have data with extremely good S/N ratio because oth- 
erwise the noise level is too high to allow the zero points caused 
by the rotation profile to show up in the Fourier transform of the 
spectral data (Simon-Diaz & Herrero |2007I >. With the rather low 
S/N ratio of our FEROS data of the faint object S 23, the pow- 
erful Fourier analysis is thus not applicable. Instead, we have to 
rely on the less accurate full-width at half-maximum (FWHM) 
method to extract the possible stellar rotation of the star, i.e., we 
fit synthetic line profiles to the observed, normalised Mgn line. 
The stellar rotation profile is described in the usual way (see, 
e.g., Gray 1976) with a limb-darkening coefficient of e = 0.6. 

The line profile of the Mgn /1448 1 line appears to be asym- 
metric (see Fig.[4j». However, this asymmetry is probably not 
real, but an artefact due to noise since only one of our spectra 
shows this asymmetry, while the other spectrum shows a sym- 
metric line profile. For a proper FWHM determination we there- 
fore rely on the blue part of the line, only. Given the low S/N 
ratio, we derive a FWHM of the observed Mgn /14481 line of 
3.5 + 0.1 A. If we assume either a pure Gaussian profile or a 
pure rotation profile, we thus obtain maximum velocity values 

Of V: 



gauss 



(140 + 5)kms and vsini =s (150 + 5)kms , respec- 
tively. However, both possibilities, i.e., either a pure Gaussian or 
a pure rotation profile, do not deliver satisfactory fits to the ob- 
served line profile (see left and middle panel of Fig.|4]i. While the 
edges are too sharp in the case of pure stellar rotation, the wings 
of the pure Gaussian are way too broad. Instead, a combination 
of both might be a good solution. But how can we estimate the 
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individual amounts of the Gaussian and the rotational compo- 
nents? 

It has been known for a long time that within the atmo- 
spheres of OBA supergiants some macroscopic line broadening 
occurs that is of no rotational origin. This broadening compo- 
nent can be described with a Gaussian profile and is usually re- 
ferred to as (macro-)turbulence (see, e.g., Howarth[2004). This 
'extra' broadening of the photospheric lines is thereby not re- 
stricted to hydrogen and helium lines, but also appears in metal 
lines. The physical understanding and interpretation of this com- 
ponent is, however, still missing. In an attempt to determine the 
contribution of this macro-turbulence to line broadening, Ryans 
et al. ( 120021 1 investigated a sample of B-type supergiants trying 
to disentangle the contributions of rotation and macro-turbulence 
with high precision. One of the key lines in their investigation 
was the Mgn 1448 1 line. Inspection of their Table 3 reveals that 
there is an obvious trend between the maximum possible tur- 
bulent velocity and the effective temperature of the star. This is 
shown in Fig. [5] where we plotted the maximum turbulent ve- 
locity found for the sample supergiants of Ryans et al. (2002 ) 
versus spectral type of the stars. While for the early-type super- 
giants maximum macro-turbulence as high as 70 km s _1 might be 
present, it rapidly decreases towards the late type stars to about 
30-35 km s _1 only. If this trend continues for the A-type super- 
giants, we might expect a maximum macro-turbulence of about 
30 km s _1 acting on the photospheric lines of S 23. 

We use this value of 30km s for the Gaussian profile and 
calculate the Mgn 1448 1 line profile by adjusting the rotational 
velocity such that the observed FWHM is reproduced. Since 
the Gaussian velocity contribution is rather small, it does not 
significantly influence the total line width, which is still purely 
rotationally dominated with about the same value of usin; ~ 
150 + 5kms _1 as without macro-turbulence. The only influ- 
ence of the Gaussian component concerns the line wings (see 
right panel of Fig.|4]i. Overall, this combination of rotation plus 
macro-turbulence gives a reasonably good fit to the observed line 
profile. 



The value of v sin i found from our fitting is high when com- 
pared to the projected rotation velocities of normal galactic B 
and A-type supergiants, for which a tendency of decreasing ro- 
tation velocity with effective temperature is reported. This can 
be seen in Fig. [6] where we plot the results for four samples: (i) 
galactic B-type supergiants from Abt et al. (2002); (ii) galactic 
A-type supergiants from Verdugo et al. ( 1999 ); (iii) SMC B-type 
supergiants from Hunter et al. ( 120081 ); and (iv) SMC A-type su- 
pergiants from Venn ( 11999I ). The SMC stars turn out to have even 
smaller velocities than their galactic counterparts, but this effect 
might arise simply from small number statistics. Compared to 
the mean rotational velocity of the galactic and SMC late B-type 
and early A-type supergiants, the value derived for S 23 is clearly 
a factor of three higher. 

However, as discussed already in the introduction, the B[e] 
phenomenon seems to be linked to rapid stellar rotation. For 
the class of the B[e] supergiants, only for three stars the value 
of v sin i could be derived so far. These are the LMC B[e] su- 
pergiants HenS93 with vsini = 65kms _1 (Gummersbach et 
al. [T995l ) and R66 with vsini SOkms 1 (Zickgraf [2006J, 
and the SMC B[e] supergiant R50 with vsini =s 150kms _1 
(Zickgraf 2000). The latter one has a projected rotation speed on 
the same order as S 23, which is now the fourth B[e] supergiant 
with confirmed high projected rotational velocity. From our de- 
rived set of stellar parameters we can further estimate the crit- 
ical velocity of the star, which turns out to be on the order of 
200km s _1 . This means, that S23 is rotating with at least 75% 
of its critical velocity. 

Interestingly, the optical spectrum of ZSW indicates that 
the Mgn 1448 1 line is less broad during the 1989 observa- 
tions. A tentative derivation of its FWHM delivers a value 
of (2.6 + 0.2) A. Neglecting any other broadening mechanism, 
this FWHM value results in a maximum rotational velocity of 
v sin i st 110 + 10km s _1 . This result is surprising, because it 
indicates that S 23 might not only have cooled, but at the same 
time also increased its projected rotational velocity by at least 
35%. 
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6. Discussion 

The investigation of our high-resolution optical spectra in com- 
bination with the photometric data of Massey (2002) resulted 
in the tentative classification of S 23 as a rapidly rotating Al lb 
star. On the other hand, our evaluation of the Mgn/Hei line ra- 
tio derived (even though with large error bars) from the spec- 
tra of ZSW, delivered a spectral classification more consistent 
with a star of roughly spectral type B9Ib. The most striking 
differences are thus an apparent temperature decrease by about 
1500 K, based on the spectral type versus effective temperature 
calibration of Evans & Howarth (2003), and shown in Fig.[3j 
and, even more curious, an increase in projected rotational ve- 
locity by about 35% (see previous section) within a time interval 
of only 11 years. In this section, we therefore investigate dif- 
ferent scenarios that might explain this peculiar behaviour. For 
this, we first consider some natural possibilities like stellar evo- 
lution (Sect. 16. it and mass ejection occulting the star in form 
of a shell phase as is often observed for the classical Be stars 
(Sect. 16.2] ). Since both attempts fail to explain the observed pecu- 
liar variations in S 23, we present and discuss in Sect. |6.3| further 
scenarios, which are more speculative. Of course, the amount 
and quality of our available data is definitely not sufficient for 
a proper analysis, which could lead to a discrimination between 
the different scenarios. Instead, we want to discuss qualitatively 
the expected observable effects for those scenarios that seem to 
be at least likely and worth further investigation. 

6.1. Stellar evolution 

The apparent cooling of the star in combination with a simulta- 
neous speed-up of the rotation velocity is unexpected and sur- 
prising. The first scenario we might think of is stellar evolution. 
However, stellar evolution theories predict that when the star ex- 
pands, i.e., cools, it can only decrease its rotation velocity. An in- 
crease in rotation speed during the post-main sequence evolution 
of a rapidly rotating star at SMC metallicity and with an initial 
mass of about 10 M can only occur during the blue loop, i.e., 
while the star contracts and heats up again (Maeder & Meynet 
120011 1. In addition, the timescales for such a star to cool from 
about 11 000 K to less than 9 500K is on the order of 5000 yr 
(Meynet, private communication). This is definitely longer than 
the one decade we found. The changes in the spectral appearance 
of S 23 can thus definitely not be explained by current theories 
of stellar evolution. 

6.2. Comparison with shell formations of classical Be stars 

A different interpretation of the apparent cooling might be given 
by the assumption of material ejection, e.g., in form of a ring or 
shell, leading to an (partly) occultation of the central star. Such 
phenomenon are known, e.g., from Be stars that can undergo 
a so-called shell phase (see, e.g., Hanuschik [1996; Porter & 
Rivinius 2003 ). The shell, ejected from the central star, thereby 
acts as an extinction source. Derivation of the stellar parame- 
ters from absorption lines created within such a shell can, there- 
fore, result in much lower apparent effective temperature val- 
ues than the real stellar one derived during a non-shell phase. 
Such a behaviour has been reported, e.g., for the Be shell star 
4 Herculis, for which an apparent drop in effective temperature 
by 3000 K during the shell phase was found by Koubsky et al. 
( 1997). However, contrary to our target S 23, but in agreement 
with expected physical behaviour, the rotation velocity derived 
from the shell absorption line decreased as well. A shell phase 



of S 23, mimicking a temperature decrease, thus has severe prob- 
lems in explaining our observations, even if the presence of some 
"shell-lines" (like the Fen and Sin lines) in our spectra hint to- 
wards the existence of a large amount of circumstellar material. 

6.3. Possible (speculative) scenarios forS23 

Since neither stellar evolution, nor a high-density mass ejec- 
tion in the form of a ring or shell can explain the peculiarities 
observed for S 23, we need to think of different mechanisms. 
In this section, we therefore present three alternative scenarios 
that might be worth further detailed investigation. These are (i) 
a high precession rate, (ii) an eclipsing binary, and (iii) surface 
He abundance and/or temperature inhomogeneities. All three are 
able to explain the peculiarities of S 23, and all three can be 
tested by future observations. But we want to stress again that 
these scenarios can only be considered as simplified models, 
given the quality of the available data. 

6.3.1 . Precession of the system 

One possible explanation for the apparent cooling with increas- 
ing rotation velocity might be a precession of the star with re- 
spect to the line of sight. Since we know that S 23 is rapidly 
rotating, the classification of the star strongly depends on the in- 
clination angle under which it is observed. 

The effects of rapid stellar rotation are manifold. Rapid rota- 
tion leads to a deformation of the stellar surface in the form of 
flattening. Thus, the polar regions become much hotter than the 
equatorial region, which is usually referred to as gravity dark- 
ening. Depending on the fraction of the rotation velocity to the 
critical velocity, the difference between polar and equatorial ef- 
fective temperature can exceed a factor of two. The same holds 
also for the other stellar parameters like the surface gravity, stel- 
lar flux, mass flux, and escape velocity. 

From the differences in observed rotation velocity, v sin i, we 
may conclude that the inclination of the system has increased 
between ZSW and our observations, which means, that the sys- 
tem was observed more edge-on in 2000 than in 1989. To quan- 
tify this, we take the values of vsini = HOkms -1 in 1989, and 
vsini = 150km s in 2000. Then, the difference in inclination 
follows to 

150 

sin 12000 = sin i 198 9 = 1-36 sin i m9 , (1) 

which means that the inclination angle has increased by almost 
40%. Assuming that the star was seen more or less edge-on (i.e., 
inclination angle of 90°) in 2000, it must thus have been seen 
under an inclination angle of about 47° in 1989. This conclusion 
is, however, only strictly valid, as long as the inclination did not 
cross the equatorial plane yet. 

A lower limit to the inclination angle in 2000 is given by 
the critical rotation of the star, u cl -i t 200 km s _1 . Assuming that 
the star is rotating at its critical speed, it delivers a minimum 
inclination angle of 48° in 2000, and of 35° in 1989. 

To get a rough idea of the influence of the rapid rotation in 
combination with the postulated precession on the observable 
mean effective temperature, we consider the following scenario 
for S 23 (which need not be correct, but which is only meant to 
show the basic idea of this toy model) : 

Since rotation at the critical limit is not very realistic, we as- 
sume in the following that the star was seen more or less edge-on 
in 2000 so that S 23 rotates with 75% of its critical speed, i.e., 
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Fig. 7. Shape of S23 rotating with 75% of its critical velocity. 
Also included are the assumed positions of the observers with 
respect to the rotation axis of the star in the two periods of optical 
observations. 
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tion see Kraus 2006). 



observed from this system seems not to be connected with the 
disk at all. 

Following Kraus (2006 ), we compute the surface structure of 
the star rotating with co = 0.75 (shown in Fig.|7]). For the adopted 
rotation speed, a deformation (i.e., flattening) of the stellar sur- 
face is obvious. Also included in this plot are the inclination an- 
gles under which the star might have been observed in 1989 and 
in 2000. While in 2000 most of the observed flux was emitted 
from the equatorial region with some minor contributions from 
higher altitudes, the flux observed in 1989 was emitted from an 
intermediate altitude, with contributions from both, the polar and 
the equatorial region. 

To obtain an idea about the expected differences in the ob- 
served stellar spectra, we further calculate the surface distribu- 
tion of the latitude dependent stellar parameters like effective 
temperature T e ff(ff), escape velocity v etic (ff), mass flux F m {8), and 
hydrogen desity «h(#) for a star rotating with u> = 0.75. These 
parameters are shown in Fig.|8]as a function of the co-latitude 6 
(i.e., 6 = at the pole), normalised to their polar values. They are 
all found to decrease by non-negligible amounts from the pole to 
the equator. 

In particular, we find from Fig.|8]the following behaviour of 
the effective temperature: 

T e «(6 = 47°) - 0.91 T eff (8 = 0°) , (2) 

r eff (0 = 47°) - 1.267^(0 = 90°). (3) 

Since for an inclination of ~ 47° the stellar emission contains 
contributions from the hotter as well as from the cooler regions, 
we might assume that the effective temperature at — 47° is 
reflecting quite well the real observed spectrum. This means, 
that we can calculate the polar as well as the equatorial value of 
the effective temperature, resulting in r e ff(# = 0°) 12 100 K, 
and T e ff(6 = 90°) - 8 700K. The equatorial value is somewhat 
lower than our effective temperature of 9 000 K derived from the 
FEROS spectrum. This is, however, no contradiction, because 
the real "observable" mean effective temperature under an edge- 
on inclination of the star will certainly be somewhat higher than 
the pure value at 90°, due to some contributions from higher 
altitudes. In addition, the temperature calibration of Evans & 
Howarth (2003 ) is certainly too coarse for a detailed tempera- 
ture assignment. 

In summary, the precession scenario might explain both, the 
apparent cooling with an increase in rotation velocity, just from 
an increase of the inclination axis of a rapidly rotating star from 
about 47° to about 90°. Whether such a strong precession will 
indeed occur, and what might be the reason for such rapid pre- 
cession movement, is, however, not known. 

To test the precession scenario, we would need to monitor the 
star over a period of several years to see whether the Hei lines 
re-appear (i.e., the star "heats up" again) with a simultaneous 
decrease of the projected rotation velocity measured from the 
Mgn lines. 



oj = v eq /v cr it = 0.75. Such an edge-on scenario is not in con- 
tradiction with our results from Sect.|4] where we found that the 
dust is most probably not located along the line of sight. In fact, 
for rapidly rotating stars the mass loss happens much stronger 
via the poles than via the equatorial region, and to date mod- 
els like, e.g., the wind-compressed disk (Bjorkman & Cassinelli 
1993 ) failed in confining the material within the equatorial plane. 
In addition, spectropolarimetric observations of Bjorkman et al. 
d 19981 ) for the galactic B[e] star HD 50138 have shown that the 
polarisation is mainly caused by scattering within an almost 
edge-on seen pure gaseous disk, while the huge amount of dust 



6.3.2. An eclipsing binary scenario 

A second possible scenario to explain the strong variations ob- 
served in S 23 might be that of an eclipsing binary system, in 
which one component is perhaps a B9 lb star and the other one 
perhaps a Al lb star. The disappearance of the Hei lines in our 
spectra might then be interpreted with the Al lb star occulting 
more or less the hotter B9Ib component, while the spectra of 
ZSW, in which the Hei lines were present and strong, indicate 
that the hotter component was dominating the spectral appear- 
ance of the system. 
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To test the hypothesis of S 23 being an eclipsing binary, it 
is certainly useful to perform a high S/N monitoring of its few 
available photospheric lines, (i.e., the Hei lines and the Mgn 
line). A further test for the eclipsing binary scenario would be 
the monitoring of its light curve. Since the hotter component is 
also expected to be the more luminous one, the system must ap- 
pear fainter when the B9 star is eclipsed than when the Al star 
is eclipsed. 

6.3.3. A rapidly rotating star with surface inhomogeneities 

A third quite interesting scenario is based on the assumption 
that S 23 might have a patchy surface abundance in He. The 
existence of such inhomogeneous or spotty surface patterns are 
known especially for some so-called chemically peculiar (CP) 
stars (see, e.g., Smith 1996; Briquet et al. 1200 II 12004; Lehmann 
et al. 120061 Krticka et al. 120071 ). These objects at or close to 
the main-sequence are supposed to be rigidly rotating, and the 
variations observed in their Hei line profiles are interpreted as 
abundance inhomogeneities in the form of large stellar spots that 
appear within and again disappear from the line of sight as the 
star revolves. In addition, these stars are found to possess longi- 
tudinal magnetic fields (e.g., Khokhlova et al. [2000; Briquet et 
al. 120071) . and there seems to be evidence that variations in the 
magnetic field might be the cause for the abundance inhomo- 
geneities. 

Furthermore, the stellar surfaces of some of the CP stars 
seem to have temperature inhomogeneities (see, e.g., Bohm- 
Vitense & van Dvk |1987l Lehmann et al. 2006), and even the 
need of a co-existence of abundance and temperature inhomo- 
geneities has recently been claimed by Lehmann et al. (120071) . 
even though the physical processes resulting in such patchy sur- 
face structures are not understood at all. 

S 23 is definitely not a main-sequence star, and to our knowl- 
edge, no such abundance or temperature inhomogeneities have 
been reported for a supergiant so far. In addition, nothing is 
known about a possible magnetic field on the surface of S23. 
Nevertheless, such a surface inhomogeneity seems to be an in- 
teresting scenario, at least to explain the variations in the Hei 
lines. Whether this inhomogeneity can also explain the different 
observed values of u sin i, needs to be investigated in more detail. 
At least the classical stellar rotation profile used in Sect.[5]to de- 
rive v sin i from the FWHM will in such a scenario certainly not 
be valid anymore. 

If S 23 possesses He spots, the variation of the Hei lines must 
be correlated to the rotation period of the star. Assuming that 
S23 is seen edge-on, the derived value of i>sin/ = 150km s -1 
delivers the equatorial rotation velocity. With a radius of roughly 
60 R e , this results in a period of about 20 days. To prove or dis- 
prove the scenario of surface inhomogeneities in either He abun- 
dance, or temperature, or even both, one would need to monitor 
the Hei and Mgn lines in the spectrum of S 23 over at least one 
rotational period. 

7. Conclusions 

We investigated the SMC star S 23, which was formerly classi- 
fied by ZSW as a B8 lb star, belonging to the class of B[e] su- 
pergiants. Our optical spectra show clear differences to the pre- 
viously published one. These differences especially concern the 
absence of all photospheric Hei absorption lines. Our spectra are 
therefore consistent more with a spectral classification of S 23 as 
Al lb supergiant. This classification is further supported by our 



analysis of the photometric UBV data. During this analysis, we 
could derive self-consistently the contributions of the interstel- 
lar extinction, and of the circumstellar extinction caused by the 
stellar wind. A circumstellar extinction caused predominantly by 
dust could be excluded. 

We further determined the projected rotation velocity of S 23 
from an analysis of its photospheric Mgn line. We found a value 
of v sin i =2 150 km s , which claims that S 23 is rotating with at 
least 75% of its critical velocity. Consequently, S 23 is the fourth 
member of the B[e] supergiant star group with confirmed high 
projected rotational velocity. 

The most puzzling result is the fact that within a period of 
only 1 1 years, between the two sets of spectroscopic observa- 
tions, the star seems to have cooled by more than 1500 K, and 
at the same time increased its rotation speed by about 35%. 
Such behaviour is definitely ruled out by any stellar evolution 
scenario. And also the scenario of a shell phase, as is often re- 
ported for classical Be stars, has severe difficulties in explaining 
all the observed peculiarities in S23. Therefore, we presented 
and discussed several other alternatives: a rapidly rotating star 
undergoing a high precession, the scenario of an eclipsing bi- 
nary, and the idea of He abundance and/or temperature inhomo- 
geneities on the stellar surface, as they have been reported from 
CP stars. All these suggestions can explain most of the curiosi- 
ties of S 23. However, the question, which of these scenarios is 
the correct one, or whether even all scenarios have to be dis- 
carded, could not be solved. Nevertheless, we suggested obser- 
vational tests that can be performed to distinguish or discard the 
different scenarios. Definitely more observations with a much 
better time resolution are needed to disentangle the mystery of 
the peculiar variations in S 23. 
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